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Abstract 



.^ • Wc analyze the approximation by radial basis functions of a hypersingular integral equa- 

tion on an open surface. In order to accommodate the homogeneous essential boundary 
condition along the surface boundary, scaled radial basis functions on an extended surface 

r^ , and Lagrangian multipliers on the extension are used. We prove that our method converges 

j^ ' quasi-optimally. Approximation results for scaled radial basis functions indicate that, for 

highly regular radial basis functions, the achieved convergence rates are close to the one 
of low-order conforming boundary element schemes. Numerical experiments confirm our 
conclusions. 
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2 '■ 1 Introduction 



This paper is about the approximation by radial basis functions of functions in Sobolev spaces 

subject to a Dirichlet boundary condition. To the best of our knowledge this is the first time that 

such functions are used and analyzed for problems with essential boundary condition. There 

5^ I arise several difficulties when imposing or analyzing trace conditions for spaces of radial basis 



functions (RBF): 

(i) Radial basis functions are selected by their center points on the domain of interest. Their 
supports are generally large and overlap with those of several other RBF. Considering 
boundary traces of RBF close to the boundary, their shapes and supports on the boundary 
vary continuously with the position of the center on the domain. Therefore, the structure 
and basis functions of trace spaces is not fixed, neither is there a fixed intrinsic basis on 
the boundary. 
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(ii) Analysis of stability and approximation properties of RBF and related operators is based 
on Fourier analysis in the so-called native space. Depending on the choice of RBF, this 
native space is a Sobolev space of high regularity. Considering domains with, e.g., Lipschitz 
boundary there is no obvious way to employ arguments from native spaces to traces. 

(iii) Traditional arguments from finite element analysis (like locality and equivalence of norms 
in finite dimensional spaces) are difficult to apply to RBF for their very nature. That 
is why the native space is of central importance. In trace spaces of RBF, however, such 
arguments (locality, equivalence of norms) are even harder to come by due to the varying 
structure of functions, cf. (i). 

In this paper we propose a mixed method employing RBF and finite elements (as Lagrangian 
multipliers) to deal with essential boundary conditions. We analyze non-conforming approxima- 
tions with scaled radial basis functions in a fractional order Sobolev space with homogeneous 
essential boundary condition. The underlying model problem is the hypersingular integral equa- 
tion on an open surface for the solution of the Laplacian in the exterior domain (with Neumann 
boundary condition). The energy space of this problem is a Sobolev space of order 1/2 on the 
surface with the condition that functions can be extended by 0. 

The analysis comprises two principal problems. One is the approximation theory for radial 
basis functions in fractional order Sobolev spaces with essential boundary condition; the other 
is the necessity of a non-conforming approach in a fractional order Sobolev space. 

Radial basis functions are well studied, mainly for the interpolation of scattered data but 
also for the approximation of partial differential equations. For some overviews see, e.g., [3 
[2 H 122 ^ 130]. In particular, Wendland studies the approximation for second order equations with 
Neumann boundary condition [29], and multiresolution properties of scaled radial basis functions 
|31j . Neumann boundary conditions allow for extending the approximation analysis to the full 
space where standard arguments apply (using the native space of the radial basis functions). To 
the best of our knowledge there exists no analysis for boundary value problems with essential 
boundary condition. In this paper we tackle this problem for a Sobolev space of order 1/2. 
Let us note that spherical radial basis functions (on the closed sphere) in this space are well 
analyzed, see |20 tll9l[25ll26j . In this paper we consider the case of a general open smooth surface 
with the particular problem of incorporating an essential boundary condition. 

This boundary condition appears in a natural way when dealing with boundary integral 
equations on open surfaces. In the case of a Neumann problem the unknown of the integral 
equation is the jump across the surface of the solution to the exterior problem [24]. Since the 
jump of this solution vanishes at the boundary of the surface the underlying energy space has to 
incorporate this condition. However, there is no well-defined trace operator in the corresponding 
Sobolev space H^f^ so that this condition appears as part of the norm; the corresponding space 
is denoted by H^''^, sometimes also referred to as Hq^ . Conforming approximations require that 
approximating functions vanish at the boundary of the surface. This causes no difficulty when 
using piecewise polynomials; the corresponding method is called the boundary element method 
(BEM). When using radial basis functions, however, a conforming and converging method is 
difficult to construct since conformity requires that the centers of the functions stay away from 



the boundary. This requirement is unreahstic because in practice centers can stay close to 
or even on the boundary. We thus propose a non-conforming approach where the essential 
boundary condition is implemented by a Lagrangian multiplier. For the BEM such procedures 
have been studied, with resulting almost quasi-optimal convergence, see [lOl [HI [H] . Here we 
use a different approach where the surface is extended to a larger surface so that supports of 
radial basis functions remain inside the extended surface. We then use a Lagrangian multiplier 
on the extended part of the surface to make the approximating functions vanish there. This 
idea is similar to a penalty or fictitious domain method. 

An overview of the rest of the paper is as follows. In the next section we recall Sobolev 
spaces, present our model problem and give an equivalent mixed formulation which will be 
used for the discretization with radial basis functions. In Section [3] we introduce scaled radial 
basis functions and the discrete scheme, and we list our theoretical results which are proved in 
subsequent sections. In Section [4] we prove the quasi-optimal convergence of our scheme, and in 
Section [5] we present an approximation theory for scaled radial basis functions. In Section [6] we 
resume our theoretical results and conclude that the resulting convergence order tends to the 
one of a standard BEM when the regularity of the radial basis functions grows. This conclusion 
is based on numerical evidence of boundedness of an additional stability term that arises in our 
analysis. Finally, in Section[7]we report on some numerical results that underline the convergence 
properties of our method. 

Throughout the paper, the symbol < will be used in the usual sense. In short, a{h, k, r) < 
b{h, k, r) when there exists a constant C > independent of discretization or scaling parameters 
h, k, r (except otherwise noted) such that a{h,k,r) < Cb{h,k,r). The double inequality a < 
6 < a is simplified to a ~ 5. 

2 Model problem and mixed formulation 

In order to introduce our model problem and its mixed formulation we need to recall the defi- 
nition of some Sobolev spaces. 

We consider standard Sobolev spaces of integer order and define fractional order spaces by 
interpolation, using the real K-method, see [2]. For a Lipschitz domain Q and < s < 1 we use 
the spaces 

H^Q) := [L2{n),H\n)]s, H'{n) := [L2(f^), ifo'(^)]. 

where the norm in Hq{Q) is the i?^(r2)-semi-norm. For orders s > 1, the spaces are defined 
by interpolation between L2{^) and a correspondingly higher integer order Sobolev space. For 
s € (0,1/2), II • ll/^sm) and || • ||j^s(r2) are equivalent norms whereas for s S (1/2,1) there holds 

H'^{i}) = Hq(^D,), the latter space consisting of functions whose traces on dO, vanish. Generally, 
H^{Q) consists of functions which are continuously extendable by zero onto a larger domain. 
For s > the spaces H~^{i}) and H~^{Q) are the dual spaces of H'^{Q) and H^{Q), respectively. 
Similarly we define Sobolev spaces on surfaces. 

In the analysis we need some more norms. In the literature different definitions of Sobolev 
norms are being used and we have to be careful to check their equivalence when combining 



different results. Apart from interpolation norms, on a domain S C IR (S = IR is allowed), 
we also need the Sobolev-Slobodeckij norm. For s = m + a with integer m > and a € (0, 1) 
we define 

„ „ /n 1,2 , ,2 V^^ 

II^IU.s := ^II^^IL,s + bL,sj 

with semi-norm 

\a\=m 



and multi-index a = (01,02)- Here, || • ||m,s is the standard Sobolev norm, as before, 

|2 
IL2(S)- 



\l,^ = E 11^"^'^ 



|a|<m 

Our model problem is as follows. For given f € H~^''^(T) find cj) € H'^''^{T) such that 

Here, F is a smooth open surface with piecewise smooth Lipschitz boundary and n is a normal 
unit vector on F. We refer to [23] for the setting of this model problem and the error analysis 
of its conforming approximation by boundary elements. 

Later we will extend F to a larger surface and use the notation W for the hypersingular 
operator on the extended surface as well. 

For ease of presentation we restrict ourselves to the case of a flat surface and assume that F is 
a polygonal Lipschitz domain in IR . Then, in particular, F satisfies an interior cone condition. 

We intend to approximate the solution (f) to ()2.ip by radial basis functions. The conformity 
condition that approximation spaces be subspaces of iJ^'^(F) requires that discrete functions 
vanish at the boundary 7 of F. This means that points for the definition of radial basis functions 
could not be freely selected since one wants to take radial basis functions of uniform radius which 
cannot be too small, as is well known. In order to be able to freely choose center points we 
extend the domain F to a larger domain F so that, for a given parameter r > there holds 
dist(7,7) > r. Here, 7 denotes the boundary of F and we assume that F is at least Lipschitz. 
Later, the parameter r will be the scaling parameter for the radial basis functions. Scaling the 
basis functions appropriately, we will be able to select center points anywhere on F and the 
discrete spaces of radial basis functions will be subspaces of H^i'^iT). 

A standard variational formulation of ()2.ip is: find cp G H^''^{T) such that 

{Wc^, ij)r = if, V')r VV' G H'/\T). (2.2) 

However, to give the setting for the discrete method we consider a non-standard formulation on 
the extended domain f: find ((^, A) G H^/^{t) x ^^^/^(F'^) such that 



ifiAW^ = V/x G #7'/'(F'=). 



(2.3) 



Here, F'^ := r\r with boundary 7U7 consisting of two connected components (T^ is an annular 
domain). Also, we use the following notation for the dual space: 

where 

with Hq-(T'^) being the space of //^(r'^)-functions whose traces on 7 vanish. In particular, any 

xp € H~' (V^) is extendable by zero to a function ip^ G H^'^{]R \ T). Below, we will use this 
notation, instead of F'^, also on extended domains F^ depending on a mesh parameter k. 
Note that H^/'^{T) = i7y^(F) and 

H-^/\T) = (^0/'(F))' = {H^/\T)y = H-^/\T). 

Also, for / defined on F, f^ denotes the extension of / by onto F (we will use this generic 
notation throughout for the extension to any domain which will be clear from the particular 
situation) . 

There obviously holds the following result. 



Proposition 2.1 Let f € H-^/'^{T) so that f £ H-'^/^{r). Then the formulations ^2^ and 
(j2.3p are equivalent and have unique solutions. There holds (f)\r = 4> and 4>\Y\f = 0, i.e. 4> = (fP , 
and A = WcfP on f \ f . 

3 Discrete method and theoretical results 



We solve (j2.3p by approximating (j) by radial basis functions and A by piecewise constant func- 
tions. To this end let $ denote a non- negative radial basis function centered around x = with 
compact support -6(0, 1) {B{y, s) is the disc {x E IR ; \x — y\ < s}) and Fourier transform 

7•($)(0 = $(6^(l + IC!')"^ eeiR2 ^3_i) 

for T > 1. The parameter r is fixed throughout. We consider scaled radial basis functions 

^r{x) := r"^<l>(x/r), r > 0, x € IR^, 

so that 

|.,,(0=^(l+r2|e|2)-\ (3.2) 

Selecting a finite set of nodes X := {xi, . . . , xn} C F we define the discrete space 
Hx,r '■= span{(l)i, . . . ,(j)]\f} where (/)i(x) := <l>j.(x — Xj), z = l,...,A''. 



Since the nodes can be near to or even on the boundary of T, the supports of the scaled radial 
basis functions are not necessarily subsets of T. We extend F to a fixed larger domain T satisfying 
supp{4>i) C r for any (pi € Hx,r, any chosen set X and any r S (0, tq] where ro > is fixed. 
We also need the mesh norm /ix,s (for S C IR ) defined by 

hx,Y: '■= supdist(x,X). 

We extend T by a strip of shape-regular, quasi-uniform quadrilateral elements T of diameter 
proportional to k as indicated in Figure 13.11 We require that the minimum diameter and length 
of the smallest edge on 7 are not smaller than k. With shape-regularity we refer to elements 
whose minimum (respectively, maximum) interior angles are bounded from below (respectively, 
from above) by a positive constant less than vr. We denote this mesh by Tk and require that 
7fc is geometrically conforming ■with 7: each element T ^Tk has either one or more entire 
edges in common with 7, or T n 7 is a vertex of T and of T. The latter case happens at most 
once for each vertex of T. 

The extended domain is denoted by T^: 

f fc = interior(f U {f; T e %})■ 

In the following, we choose the mesh size k accordingly to the scaling parameter r: k > r and 
k is small enough so that F^ C F and no element touches two vertices of F. The assumption 
k > r guarantees that the supports of the scaled radial basis functions are within F^. 

We introduce the notation F^ for the strip F^ \ F and 7^ := dV^ for the boundary of F^. 

Let us collect the assumptions we have made: 

(Al) The scaling parameter r and mesh parameter k are bounded, r S (0, rg], r < fc < feo- 
7fc is a quasi-uniform mesh of shape-regular quadrilaterals of diameter proportional to k 
with minimum diameter not less than k and with edges on 7 of length not less than k. 
Moreover, Tk is geometrically conforming with 7 and no element touches two vertices of 
F. 

Now, for the approximation of the Lagrangian multiplier A we take piecewise constant func- 
tions, 

Mfc := {/i G L2{Tl); /i|r = const VT G Tk}. 

Using these discrete spaces, the boundary element scheme with radial basis functions and La- 
grangian multiplier for the approximate solution of (j2.3p is: find {(j)N,^k) £ Hx,r x Mk such 
that 

{fi,MTi =0 y^ieMk. ^ ■ ^ 

Defining the subspace 

Vx,r ■■= {^ e Hxy, ii^, i'Wt =0 V^i G Mfc}, 
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Figure 3.1: Domain F extended by strip F^ with mesh 7fc. 



and assuming an inf-sup condition for the bilinear form (•, •)rc on M^ x Hx,r, (|3.3p is equivalent 
to: find (j)N € Vx,r such that 

Apart from assumption (Al) we will need some more properties of X for our analysis. We 
must be able to fix a constant function on any given element by testing with a radial basis 
function. 

To make this precise, let C C 7fc denote the set of elements T such that the closure T 
intersects 7 only at a vertex of F. For instance, in the example of Figure [3?T] there are five such 
elements, at five of the six convex vertices. Further we denote £ := Tk\ C. Then we assume: 

(A2) Vr G ^ 3xieX : supp(0i) n F£ C T 

(A3) VT G C : T has exactly two neighbors Ti, Ta G £, 

3xieX: supp((/.i) n T / 0, supp(0i) n T = VT G Tfc \ 1^, Ti,T2} 

For any T ^ Tk there may be more than one Xj G X satisfying either (A2) or (A3). We will 
denote one of these points Xi by Xj(t) . 

Moreover, we assume that there is substantial overlap between Hx,t ctnd M^: 

(A4) 3k > : meas(supp(0j(2')) n T) > k meas(supp((^j(r))) = Kvrr^ VT G Tk 

Remark 3.1 Note that the discrete scheme is defined on the domain Tk which depends on k. 
But this causes no difficulty since, under the assumptions made, all the domains Tk can be 
extended to the fixed domain T and the spaces involved allow for extension by zero of their 
elements to T. More precisely, there holds 

Hx,r C H'^Htk) C H'/\t), Mk C H~'/\Tl) C H~'/\T'^) 



where the latter inclusions are to be understood as the uniformly continuous injections of the 
respective extension by zero. 

Theoretical results. Based on assumptions (A1)-(A4) we prove the foUowing results. 

• The discrete scheme (j3.3p converges quasi-optimally; see Theorem 14. li 

• The error of the best approximation by radial basis functions in the constrained space Vx,r 
(mean value zero on elements of the extension) can be bounded by the error of the best 
approximation in the unconstrained space Hx r on F plus a stability term; see Theorem l4.3i 



We prove an error estimate for the best approximation by scaled radial basis functions 
(Theorem 15. 4p which appears to be sharp according to our numerical results in Section [71 



We are unable to show an appropriate bound for the stability term |^|/f-t(rc) in Theorem 14.31 
Based on our numerical tests, however, we conjecture that this term is appropriately bounded 
so that we conclude for the choice A; ~ r ~ h^ p the overall error estimate 

U-<pN\\m/2(^r) < C h%-^'^-'^'^^-"^^ e > 0,t G (1/2,1]. (3.4) 

The constant C would depend on t, e and r but not on X and hx,r under the assumptions 
made. We refer to Section [6l in particular (16. 6p . for more details. 

Remark 3.2 In our case of an open surface T, the solution (p of (12. ID has strong singularities 
along 7. This limits the convergence order of approximation schemes. Measuring regularity 
in standard Sobolev spaces, the h-version of the standard boundary element method with quasi- 
uniform meshes (and mesh size h) converges like 

U-Mm/2^r)^h'-'^'MH^[r), (3-5) 

cf. |27j. Since (p £ H^i^) for o,ny t < 1 but (p ^^(r) in general there is an upper limit 1/2 
for the convergence order in h. An optimal error estimate making use of the type of appearing 
singularities is 

again for quasi-uniform meshes, see [3]. In our case of radial basis functions the mesh size h 
corresponds to the mesh norm hx,r o,nd for quasi-uniform distribution of nodes this parameter is 
equivalent to h for quasi-uniform meshes. The estimate (j3.4p exactly reflects the error estimate 
(13. Sp for large r. Selecting sufficiently smooth radial basis functions, which corresponds to large 
T, one gets as close as wanted to the convergence order 1/2. The use of lower regularity of radial 
basis functions results in a lower convergence order. Our numerical experiments reported below 
confirm the predicted influence of t. 



4 Quasi-optimal convergence 

We prove the quasi-optimal convergence of (j3.3|) for the approximation of 0. Later, in Section [5] 
we study the approximation problems for (j) and A to derive convergence orders. We do not 
bound the Galerkin error for the Lagrangian multiplier A since, on the one hand, proving a 
discrete inf-sup condition for the bilinear forms (•, •)r'= is an open problem and, on the other 
hand, the function A is of no physical interest. We will therefore analyze (j3.3p without using a 
discrete inf-sup condition. 

Theorem 4.1 Let the assumptions (Al)-(A3) he satisfied. Then there exists a unique solution 
to (|3.3p and there holds the quasi- optimal error estimate 

\\clP - 07vbi/2(f^) < ^mf ^^ 11/ - nm/Ht,) + Jf,^ l|A - ^^Wn-^i^^rir 

Proof. The existence and uniqueness of {(jiN-, ^k) G Hx,r x ^k follows from the Babuska-Brezzi 
theory. Specifically, uniformly in k, 



(i) (VF-,-)p is bounded: 

{Wv,ij)f.^ = (l^i;°,V'°)f ^ ll^°lli/i/2(f)llV'°llHi/2(f) - Mm/HrJ\'^\\m/Hr^) (4-1) 
for any i;, ^ G #^/^(ffc) 
(ii) {W-, •)p is elliptic: 

(T^V',V')f, = (^^<,Of ^ IIV'°lll,i/2(f) ^ IIV'll|i/2(p^) vv e Hy\fk) (4.2) 

(iii) (•,-)r= is bounded: 

{fi,i;)ri < ||/i||^-i/2(^.)||V'||^i/2(p.) < ||^.||^-i/2(p.)||^||^i/2(f,) (4.3) 

for any fi G #7^/^(r|) and ip G ^^/^(ffe) 
(iv) The linear form defined by f^ is bounded, i.e. 

Here, we used several uniform norm-equivalences, e.g. ||'y||j^i/2/p ^ — ||'y''||j^l/2/p^ since H^''^{Tk) 
(resp., iJ^'^(r)) is defined by interpolation between L2(rfc) ^ind HQ{Tk) (resp, between L2{r) 
and Hq{T)), and 

Mn^oi^k) = \^\HHn) = l^°lHi{f) = ll^°llHi(f) ^^ ^ -?^o(ffc). 

Rather than proving a discrete inf-sup condition for the bilinear form (•,-)r= we only show 
injectivity. More precisely, 



(v) (•, ■)ri is discrete injective: 

[fieMk-. (/x, V)r = V^ G Hx,r ) ^ ^ = (4.4) 



To see this we proceed as follows. Let fj, = YlreT- '^tXt (with characteristic function xt 
on element T) satisfy (|4.4p . For an element T ^ S, i.e. T has at least an entire edge in 
common with 7, there exists, due to assumption (A2), a basis function (/>j G Hx.r whose 
support overlaps only with the element T. Then {fi,(j)i)Y'c = CT{xT,'Pi)r'' = so that 
ct = 0. Therefore, ^ vanishes on all those elements. Elements not having an entire edge 
in common with 7 can only be at convex vertices of F and are isolated. By assumption 
(A3) we can again choose a basis function for each of those elements T, this time with the 
only condition that there is some overlap between T and the support of the corresponding 
basis functions. Since we already know that ^ vanishes on the neighboring elements, the 
argument from before implies that /i vanishes also on the remaining vertex elements. This 
proves (j4.4p . 

The Babuska-Brezzi theory then implies that there exists a unique solution {(J)n^ ^k) of (|3.3p . 

To prove the quasi-optimal convergence we first derive a Strang-type error estimate. Follow- 
ing the standard procedure, we use the triangle inequality and the uniform ellipticity (14. 2p to 
conclude that for any -0 G Vx,r there holds 

ll-^" - <?^Af lli/l/2(ffe) ^ ll'^" - V'llHl/2(f,,) + Un - V'llHl/2(ffc) 
<IUO ;„ ^ {W{cPx-i^),^)t, 

^ll/-^llffi/2rf„^+ sup 



r\^ 



H^'HTk) 



vGVx,r\{o} \m\H^/-^{fu) 



<ll/-^bl/2(f^) + 



rv-' 



(W(/-V'),V')f,. (VF(/-07v),'/')f, 



sup 7 — 7 + sup *" 



Then, using the uniform boundedness ()4.ip . we find 

This is the Strang-type error estimate for a non-conforming approximation. For a conforming 
method the latter term above vanishes. 

It remains to bound the non-conformity term. Combining (|2.3p and (|3.3p one finds that 
there holds (note that ^ = / by Proposition 12. ip 



{W{(tP - (I)n), ^)y^ = (A - Afc, ip)Ti = {X- fi, ip)ri V99 G Vx,r, V^u G M^. 
Combination with (j4.5p and application of (|4.3p yields 

11/ - 0^llHV2(f,) < ^inf ^^ 11/ - ^llHV2(f,) + ^^ l|A - Hl^-i/2(r.). 
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This finishes the proof. D 

In the next step we analyze the approximation error in the constrained space. 

Lemma 4.2 Let the assumptions (A1)-(A4) be satisfied and let s,t € (1/2, 1] with s < t. Then 
for any (p G -f^*(r) there holds 

inf ||/-V||j^.(f,-) < C(s,0 inf ( ||(/) - V'lU^cr) + (^(s,*, ^,'^)lV'lH*(rn 

with 

6{s, t, k, r) := A;i+V-(i+^) + A;*"^ (l + k^+'r-^^+'^) . 

Proof. Let ipi £ Hx,r be a ininimizer of 

\\4>-i^\\H''(r) +S{s,t,k,r)\tl;\Ht(^ri) 
among the elements ip € Hx,r- We will construct a function 1^2 £ Hx,r such that 

(/i,^i-V'2)r^ =0 yfieMk, (4.6) 

that is ■0 := V'l ~ "02 £ ^x.r, and such that 

11'^° - V'lli/.(f ,) < C{s, t){\\(l) - V'l lli?=(r) + <5(s, t, k, r)\i;i \miri)) ■ (4.7) 

To this end recall the notation of C for elements that touch 7 only at a vertex and of £ for 
elements touching 7 with at least an edge. Let us denote hj £0 C £ the set of elements which 
do not touch any element of C. 

By assumption, the radial basis functions are non- negative. We consider the normalized 
functions (p* := 4'i/\\4'i\\Li(T) for i = i{T) so that ||(Ai'||Li(T) = 1- This normalization is well 
defined since there is substantial overlap between supp(i;^j(Qn)) and T by assumption (A4). We 
make the ansatz 

For T £ £q we define 

CT = / V'l 

Jt 
so that (1, V'l ~ V'2)t = due to the chosen normalization. 

If C 7^ then there remain elements in £ which are associated with vertices of T. Let us pick 
one vertex where this happens, i.e. there is an element T2 € C touching 7 at this vertex and 
there are two neighboring elements TijTs S £\£q. For illustration see Figure I^TTl 

For this vertex we define 

CTi = / V'l - / (l)*im) / V'l, CT2 = / V'l, CT3 = / V'l - / 4>*im) / ^^■ 

JTi Jn JT2 JT2 JT3 JTa Jt2 
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Figure 4.1: Corner elements and supports of associated radial basis functions. 



Repeating this construction for all elements T S C we obtain ■02 G Hx,r satisfying (j4.6 
It remains to verify (j4.7p . Since s G (1/2, 1] there holds by [3l Lemma 3.1] 



P- - 



)-|l^°-*^(f.) 



1^- 



<2||/-V||?,r + 2||/ 






Here, we used the equivalence of the interpolation norm || • Hj^s/p ■> and the Slobodeckij norm 
II ■ ll^p with constants independent of k due to meas(r) < meas(rfc) < 2meas(r). However, 
they may depend on s and may be unbounded when s — > 1/2. By noting that cj)^ = (p on F and 
(p^ = on Vf,, and by using [H Lemma 3.5] we deduce 



l^nrfc 



T€% 



,)<II'^-V'll?,r+E(^"''II^IIL(T) + IV'I',T)- 



Using the triangle inequality with ip = ipi — tp2 and the equivalence of norms again, this time 
on r, we obtain 






^ \\(t> - ^illl/-(r) + \\M\'hs{v) + ^ (k 



T£% 



-2s I 



\L2{T) 



+ 



\It 



Now, the support of 1^2 is confined to a neighborhood of the boundary of F. Therefore, by a 
Poincare-Friedrichs inequality its norms can be replaced by the semi-norm, giving 



|2_ . < 



Mh(^r) + l^2|^.(r) +Y.[^ ''ll^llL(T) + IV'I'.T 

TeTk 



We finish the proof of the lemma by showing that 

lV'2b«(r)<A;^+V-(i+^)|V'i|^.(r^) 



(4.8) 
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and 

E (^-'iiv'IIL(t) + i^i^,t) <^^^*-^)(i+^^+v-(i+*))^i^,i^.(p.). (4.9) 

TeTk 
Proof of (|4.8p . By a coloring argument and the Cauchy-Schwarz inequality we start bounding 

IV'2|^s(r) < E 4l'/'*(T)ll^=(supp(0* ))• (4-10) 

TeTfc 

Here we used that only a fixed number (independent of all relevant parameters) of appearing 
radial basis functions overlap. By the scaling property of the ff**-semi-norm (see, e.g., |13j ) there 
holds 

|i |2 2~2s\i |2 2— 2s —4 

and by the assumption of substantial overlap (A4) one finds 

II4(T)IIli(T) - 1- 



This proves 

l<^*(T)l//-(supp(</,*(y))) = ll'/^i(T)llLi(T)l'^»(T)l//-(supp(</,*(y))) 



*,„J2_, _ ,. = IU.,^J|72^._,U..^J2_ , _ ^ _^-2(l+s)_ (4_^;^) 



Now, for T (^ Sq, again using scaling properties, transforming to a reference element T, de- 
noting the transformed function by adding the symbol " " " , and applying a Poincare-Friedrichs 
inequality, we obtain 

4 = ij^i^if ^ kHJj^f < k'\iPi\%^^^ c. k'+^^\^,\l,^^y (4.12) 

For a corner element T2 € C with neighboring elements TijT^ (^ £ we obtain 

(4.13) 





4,<t'«'l*.l?,,7.,, 


as before and 




4,<2{l*>)\2{ 


i/-'r^>i,*-^y '-''*'>' 



In the last step we used that Jj, (f)*,j. -, ^ 1 by the quasi-uniformity of the mesh, the normalization 
/t ^*i{T ) ~ ^ ^"^^ ^^^ substantial overlap of supp(0*/;j, J with T2. Accordingly one bounds 

and repeats this procedure for all edges where necessary. Combining (j4.1ip - ()4.15p and recalling 
(|iTO]) we obtain (1^8]) . 
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Proof of (j4.9p . We use scaling arguments and a Poincare-Friedrichs inequality as before. By 
the integral-mean zero condition (or using that ip satisfies a homogeneous boundary condition) 
and scaling properties we can bound 

and 

so that 

E (^"'IV'IIL(T) + IV'I^.t) < E k^'-'-'^mn^iT) < k'^'-'WH^^r^,y 

Analogously to (14.80 we can bound 

lV'2|Hnr^.)<A:^+V-(i+*)|Vi|H*(r^,)- 

Therefore, the representation ip = ipi — ip2 and the triangle inequality, together with the previous 
estimate, yield 



This finishes the proof of (j4.9p . D 

We can now state the main result of this section, the quasi-optimal convergence of our scheme 
(|3.3p in the unconstrained space Hx r x Mk ■ 



Theorem 4.3 Let the assumptions (Al)-(A4) be satisfied and let t € (1/2,1]. There exists a 
unique solution (p to (j3.3p . If (j) £ H^(X) then there holds for e > the error estimate 

ll/-'/'7vllHi/2(ffc)<C'(e) in^ (||(/)-V'||j:^i/2+.(r) +(5(e, t, /c,r)|^|^t(pc)j + inf ||A-^||^-i/2(p,^ 

where C{e) > depends only on e and with 

6{e, t, k, r) := fc^+V-^/^"^ + k'-'/^-' (l + k'+'r-^'+''>) . 

Proof. The statement is a combination of Theorem 14.11 and Lemma 14.2^ setting s = 1/2 -|- e in 
the latter one. D 

In the next section an error estimate for the best approximation by scaled radial basis 
functions will be derived. 
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5 Approximation by scaled radial basis functions 

Standard approximation theory of radial basis functions centers around the native space defined 
by 

Afq> := \v e L2(M'^); / ) d^ < oo j with norm \\v\\j^, := \\^~^^'^v\\. ,-rn2. 

which, for ^ satisfying (j3.ip . is identical to the Sobolev space H'^[]R ) with equivalent norms 
when defining 

IIHI^j(]R^):=ll(l + |-Pr/'^(-)llL.(]R^) («elR). 

It is well known that the interpolation operator Ix defined by 

Ixv e Hx =span{^{- - Xi); i = 1,... ,N}; Ixvixi) = v{xi), i=l,...,N 

for V £ A/i, satisfies 

\\v - IxvWm^ = min \\v - iI^Wj^^, , 

i)&Hx 

see [8]. 

Using the scaled radial basis function ^j. the native space 7V$^ is still identical to //'^(IR ) 
but its norm is uniformly equivalent to the norm 

II^IIh^,(]R2)-II(1+^'|-|')'^'*(-)IIl2(]R2) (^>0) for s = r, 



cf. ([32]). 

For the analysis, we also need the scaled versions of the norms defined in Section [2l For 
s>0, s = r7T, + (T with integer ?tt, > and a G (0, 1), the scaled Sobolev-Slobodeckij norm is 
defined by 

1 /9 

\\v\\s,r,^ ■■= (\\v\tr,j: + r^"\v\ls) (^ > 0) 



where S is a domain in IR and 

l|2 _ 'V^ 2|a|||na \\2 



|2 ._ V^ 2|a|||r,a ||2 



\a\<m 

with multi-index a = (ai,Q2)- We also define the scaled interpolation spaces 

H^{J:) = [//r(S), //r+^(S)]a, is = m + a as before) 

with norm denoted by |H|h^(e)- By the equivalence of the semi- norms l^l^jp^a and |lH'^i'(-)lli m^) 
(cf. |16|. Lemma 3.15]) and by repeating the arguments of Theorem B.7 in J16j we deduce that 

Ms,rM' - II^IIh^,,{]R^) - MmiM') ^^ e H'OR^), r > 0, (5.1) 

where the constants are independent of r. In particular, when s = t there holds 

IM!r,r,]R2 - MlHrQR^) - Iblk*. Vt; € H^OR^), r > 0. (5.2) 

Furthermore, we need the following result. 
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Lemma 5.1 For any s > there exists a bounded extension operator 

E: H^{r) ^ H^{M^). (5.3) 

As a consequence, there holds 

\\v\\s,rT < \\v\\H?ir) V« G H%T). (5.4) 

In both cases the boundedness is uniform for r > bounded from above. 

Proof. By Stein (see [231 Section 3, Chapter VI]) there is a bounded extension operator 
E : H''^{r) -^ H"^(\R^) defined for ah non- negative intergers m. It follows that E : H^{T) -> 
H^(JR ). Indeed, for any integer m there hold 

|Q!|<m |a|<7« 

~ Z^ '^ ll^lli/l-l(r) - Z^ ^ 2^ ^ 11^ ^llL2(r) 

|«|<m |a|<7?i |/3|^|c* 

^ Z^ Z^ ^ W^ ^\\L2ir) ^ Flln^cr)' 

|o|<m|/3|<|a| 

where in the penultimate step we used the fact that r is bounded above. By interpolation we 
obtain the boundedness of E : H^{T) -^ H^(\R^) for ah s > 0, i.e., (fOjl . 

To prove (|5.4p we note that Ev = u on F to obtain for any v € H'^{T) with s = m + a 



Hlr. 



ta|<m |o|<m 

\a\<m 

By using ([53]) and (fOj) we deduce ([531) • □ 

In the following we recall and adapt techniques from [18] to bound the approximation er- 
ror appearing in Theorem 14.31 We will need the following result; see |18t Theorem 2.12], [TJ 
Corollary 4.1]. 

Proposition 5.2 Let k be a positive integer and a € (0, 1]. Then there exists Hq > such that 
for any X with /ix,r ^ ho ^iT'd for any I = 0, . . . ,k there holds 

\v\i,r < C{k,a)h''+^-^\v\k+a,r ^v G i/'=+'"(r) with v\x = 0. 

In the following, [rj denotes the largest integer smaller than or equal to r. 
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Lemma 5.3 Suppose that assumptions (Al) and (j3.ip are satisfied. Then for s G [0, [rj] there 
holds 

\\v-Ixv\\H?(r)<Cis,T)i—^j \\v\\H-{r) Vre(0,ro]. 

Proof. Let 

E: H;{T) -^ H^OR^) 

be a uniformly (in r) bounded extension operator, cf. Lemma l5.ll and let v € H^{T). 

Using that Ev = t; on F and thus IxEv = Ixv = EIxv on T (since X C T) one finds 
that Ev — IxEv is an extension of u — Ixv. Therefore, the property that Ix is an orthogonal 
projection in A/^^ yields, noting (j5.2p . 



\\v - IxvWh^^t) ^ ll^fc^ - IxEkv\\fjr(jR-2) ^ \\EkV - IxEkv\\js[^^ 

< \\Ekv\\M<,, ^ ll^^''^ll/f,T(IR') - Mn-ir)- (5-5) 

For integer £ < [rJ , Proposition 15. 2| estimate (|5.4p and stability (j5.5p yield 
i\ 1 \ ^ e.'LT-e.w T II ^ ( ^x,t\ II J II 



^ ,'/ix,rV II r II / /^^x,rV II II 



This proves the assertion for integer s. For non-integer s < r we interpolate between i/j. -'(F) 

and Hr (F), noting that if s = [sj +0" with a G (0, 1) be such that s < [rJ then [sj + 1 < [t\ . 
D 

Theorem 5.4 Xei assumptions (Al) and (j3.ip 6e satisfied. For t G (0, \t\] with t + 1/2 being 

non-integer, let r ~ h-^-p . Then for any v G H^{T) there exists ip G Hx,r such that for 
< s < t there holds 

(t-s)il-t/T) . 



\v -i^Wn^'ir) < Ch 



V 



Here, the constant C is independent of v and /ix,r ^ut may depend on s, t and r. 

Proof. We follow the proof of jl81 Theorem 3.8]. Let v G -f^*(F) be given. According to 
[171 Proposition 3.6] for any o" > 0, there exists a band-limited function g^^ G B^^ := {v G 
L2(]R^); supp({)) C B{0,a)} such that (noting that || • ||^i,(]f^2 ~ || • ll^s^j^aJ 

\\v" - 9a\\Hs(R^) < Cis,t)(T''~^\\v^\\Ht(^Si2y (5.6) 
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We then define '0 '■= Ixda and obtain, together with Lemma 15.31 

\\v - Ixga\\Ht{r) ;$ W'v - 5'(t|Ih*(]r2) + \\9o- - Ixga\\Ht{r) 

^ ll^° - 9a\\Ht{We) + r~*'\\gcT - Ixga\\Ht{r) 



- Mh^{T) + r'^h^'^WgAn-ir)- (5-7) 



Here, and in the rest of the proof, h denotes hx,r- Also we used that 

when t + 1/2 is not an integer, cf. J15j . 

Since g^ is a band-limited function there holds 



< (1 + ,V)-* /■ |^.(0P(1 + r^lep)* dC^{l + r'aT%a\\'.y 

Again using (j5.6p . we obtain 

II5'ctIIh*{]r2) ^ II5'(7||/^*(]r2) < ||u - gu\\Ht{we) + H^ IIh*(ir2) ~ 11^ IIh*(]r2) - ll^llift(r)' 

so that with the previous estimate. 

Combination with ()5.7p yields 

11^ - /x5<x||H*(r) < (l + r--/i^'*(l + rV2)(-*)/2) ||^||^,(^^. (5.8) 

By the same arguments and using the same construction we bound 

11"" - ixga\\L2{T) < 11^ - gAi-ii^^) ^ "^'^ ~ ixgAL2{T) 

= (^-t + ^-/,-(i + rV2)(--*)/2)||^||^,(^^. (5.9) 

Now choosing o" = 1/r in (j5.8p and (j5.9p . and taking r = h"' with a = 1 — t/r, we obtain 
Ik - -^A'9<7|k*(r) < (l + r~^/i^"*) |l?^|lj^t(r) = 2 |lf |l^t(r) 
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and 

Interpolation gives 

11^ - Ixga\\H'^(r) ^ \\v - Ixga\\Ht{r)\\'" ~ -^^^''^IIl^ct) 

This proves the theorem. D 

6 Conclusions 

Before presenting numerical experiments let us draw some conclusions. 

Based on Assumptions (Al)-(A4) we have proved the quasi-optimal convergence of the 
discrete scheme (j3.3|) (cf . Theorem 14. 1|) . The best approximation of the Lagrangian multiplier is 
taken in the natural space of order —1/2 on the domain of definition of the Lagrangian multiplier 
(outside r). The best approximation of the sought solution (p is measured in the natural space 
of order 1/2, on an extended domain and taken among scaled radial basis functions of the 
constrained space, i.e. with piecewise mean value zero on the extension. With Theorem 14.31 we 
managed to replace the latter term (best approximation of (j)) with the best approximation in 
the unconstrained space on the original domain T. The price to pay is an additional stability 
term that measures the approximant on the extension, where the unknown solution vanishes. 
We were able to bound the best approximation error on the original domain by a term that 
shows a convergence order that is close to the one of a standard boundary element method when 
the regularity r of the native space becomes large, cf. Theorem 15.41 

We were unable to show an appropriate bound for the stability term |'0|H*(r'=) in Theorem l4.3[ 
The natural tool to estimate |V'|_H'«(r'=) = l-^x5o-|_H't(r'^) is switching to the norm of the native 
space. (Recall that §„ is a band-limited function that approximates the solution of the integral 
equation.) This switch produces a factor of r""^ for the L2(r)-norm of g^ which we cannot 
control efficiently as we can with the other higher order term \ga\H-^{r)- 

Our numerical results for the choice k c^ r c^ hy- p and t = 1 (cf. Figures l7.2H7.4p in- 
dicate that |V'|_ff«(r'=) asymptotically behaves, for the discrete solution (f)^-, exactly as the best 
approximation error. Note that for the chosen parameters, 

5(e,t,fc,r)|V'U.(r|) ^ k'~^'^"'\i^\miTi) ^ 4:r^'-^^^'"*^"VlH*(r^,), (6.1) 

so that it is enough to have boundedness of \iI^\h*(t'') to obtain the optimal approximation order, 
cf. the final estimate (j6.6p below. Based on the assumption that the minimizer ^ € Hx,r of 

||(/) — ^||j:^i/2+e(r) has bounded semi-norm |V'|_f/*(r'=)) i-6-, 

inf (110- V||i:/i/2+.(r) + '5(e,t,A:,r)|'(/'|//*(rn) - , inf U - M\mn+-(VY (6-2) 
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let us deduce a final error estimate. Under assumption (|6.2p . Theorem 14.31 gives 

11/ - <A7v||^i/2(f,) < CCe) , inf U - V'llHi/2+.(r) + inf ||A - /i|U-i/2 .. (6.3) 



By Theorem 15.41 we bound 



(t-l/2-e)il~t/T) , 



, inf U-nH^/.+.^r)<hTr' " ' 'MH^ivy (6-4) 

There holds A = (VF/)|r-, cf. Proposition EH Now, / G ^*(f) so that A|r- S -H'*~Hr|) by 
the mapping properties of W, cf. [16]. A standard approximation result (see [5l Lemma 2.3] for 
the p-result on an element; this immediately generalizes to the present case) yields 

Mk "^ " ^W^^W - k'-'/'\\X\\H^-^^riy (6.5) 

Combining ()6.4p and (|6.5p with k c^ h - , and using the estimate (j6.3p . this leads to 

-^1^ k 

\\4>-4>N\\m/2^r) < ll'/'°-</'iv|lHi/2(f,) ^C4:r^'"'^^'"*^"^ e > 0,t G (1/2, 1]. (6.6) 

The constant C would depend on t, e and r but not on X and hx,r under the assumptions 
made. 

7 Numerical results 

We consider the model problem (j2.ip with F = (0, 1) x (0, 1) and / = 1. The nodes of X 
are distributed uniformly on f where F is extended by a strip of uniform squares (with side 
length k) as in Figure 17.11 The support of the scaled radial basis functions is indicated in one 
case. The setup is selected so that assumptions (Al)-(A4) are satisfied and with mesh norm 
h = h^ J. c::i r c::^ k, according to (|3.4p . We use scaled radial basis functions with the radial basis 
functions defined in |28] (for the case d = 2 there which corresponds to functions in IR ). These 
functions satisfy relation (|3.ip with r = 3/2 + m where m is the parameter of the corresponding 
regularity C^™. The degree of the corresponding univariate function is 2 + 3m. 

We calculate the errors in an approximating (and heuristic) manner. For a conforming 
method, by making use of the symmetry of the hypersingular operator, one obtains 

11/ - '^^ll|i/2(f) ^ (w^(/ - M, (/ - .^^))f 

= (T^/,/)f - {Wcl)N,<pN)r = {W(l),cl))r - {W<PN,Mr- 

The last term is available through the stiffness matrix of the problem and the term {Wcp, (/>)r 
can be approximated by extrapolation, cf. [9]. In our case of a non-conforming (or mixed) 
approximation, this calculation has a perturbation which is due to the term A = (TV (/>'') |r<= 7^ 0. 
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Figure 7.1: Uniformly distributed nodes X on F, and uniform mesh generating F^ and extending 
Ftof. 

Since we do not know the exact solution ^ so that A cannot be calculated we approximate the 
relative error in the energy norm by the expression 



lex 



(VF(^Ar,(/>Af)p 



1/2 



(7.1) 



Here, ||</>||ex denotes the extrapolated value substituting 1^^^(^)y- 

For different values of r, we present the approximated errors on double logarithmic scales 
along with the expected convergence rates (upper limit) according to (13. 4p . For comparability 
we use the same scales in all the figures. Table 17.11 lists the values of r with corresponding 
figure number, data (regularity C^*" and polynomial degrees as mentioned before) and the limit 
1/2(1 — l/r) for the expected convergence rates. 



r figure no. 



3/2 
5/2 
7/2 



E2] 



El 



regularity polynomial degree exp. conv. rate 2(1 — 7) 






1/6 
3/10 
5/14 



Table 7.1: Values of r in the numerical experiments with corresponding data, expected conver- 
gence rates and figures. 



Figures I7.2H7.4I confirm quite precisely the predicted convergence. The lines indicated by 
"error" give the approximate relative errors calculated by (17. ip whereas the lines "stab term" 
give the values of the stability term (5(e, i, fc, r)|'i/'|H*(r'= for e = 0, t = 1 (the limit of the regularity) 
and V the calculated RBF approximation, cf. (16. 1|) . The lines "expected" plot the expected 
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convergence rates as listed in Table 17.11 They correspond to our conclusion ()6.6p for e = and 
t = 1. All results indicate that the errors have the predicted convergence rates and that the 
stability term (j6.ip fulfills our conjecture (j6.2p . 

We have implemented the method by numerical integration with an overkill of number of 
integration nodes. We used transformation to polar coordinates so that the singularity from the 
fundamental solution cancels in the diagonal entries of the stiffness matrix. Nevertheless, note 
that the polynomial degrees for larger values of r are large (8 in the case t" = |) which makes 
their implementation a non-trivial task. 

In the case t = | there is a large pre-asymptotic range (Figure [73]) • Note also that for t = | 
we were not able to calculate the stability term for the whole range of unknowns (about 30,000). 
In this case the radial basis functions are only continuous and the numerical calculation of the 
if ^-semi-norm becomes unstable. 
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Figure 7.2: Relative error and theoretical convergence rate for r = 1.5. 
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Figure 7.3: Relative error and theoretical convergence rate for r = 2.5. 
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Figure 7.4: Relative error and theoretical convergence rate for r = 3.5. 
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